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ABSTRACT 

We consider nucleosynthesis in outflows originating from the inner regions 
of viscous accretion disks formed after the collapse of a rotating massive star. 
We show that wind-like outflows driven by viscous and neutrino heating can ef- 
ficiently synthesize Fe-group elements moving at near-relativistic velocities. The 
mass of ^^Ni synthesized and the asymptotic velocities attained in our calculations 
are in accord with those inferred from observations of SN1998bw and SN2003dh. 
These steady wind-like outflows are generally proton rich, characterized by only 
modest entropies, and consequently synthesize essentially nothing heavier than 
the Fe-group elements. We also discuss bubble-like outflows resulting from rapid 
energy deposition in localized regions near or in the accretion disk. These inter- 
mittent ejecta emerge with low electron fraction and are a promising site for the 
synthesis of the A = 130 r-process peak elements. 

Subject headings: gamma rays: bursts — nucleosynthesis — accretion disks 

1. Introduction 

In this paper we examine the production of nuclei in matter escaping the innermost 
regions of collapsars. Collapsars occur when the usual neutrino-powered supernova shock 
fails to expel the mantle of a rotating massive star whose core has collapsed (Woosley 1993; 
MacFadyen & Woosley 1999). The inner parts of the star collapse to form a disk accreting 
rapidly onto a central black hole. Interesting elements, and in particular ^^Ni and A = 130 
peak r-process elements, may be synthesized in outflows from the inner regions of this disk. 
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Understanding how coUapsars make ^^Ni, whose decay fuels optical light curves of SNe, 
is central to connecting the deaths of massive stars with Gamma Ray Bursts (GRBs). There 
are now observations of SN-like light curves for some five or six GRBs (Price et al. 2003). 
For two SNe/GRBs (SN 1998bw-Galama (1998) and SN2003dh-Stanek et al. (2003); Hjorth 
ct al. (2003)) there are detailed estimates of Ni ejecta mass and velocity (Patat ct al. (2001); 
Iwamoto et al. (1998); Woosley, Eastman, & Schmidt (1999) for SN1998bw and Hjorth et 
al. (2003); Woosley & Heger (2003) for SN2003dh). 

Though it is generally agreed that collapsars are promising sources for the observed Ni, 
details remain uncertain. There arc several possibilities. As in "ordinary" SNc, Ni may be 
synthesized explosively as a strong shock traverses the stellar mantle. Parametrized piston- 
driven simulations of the explosion of massive stars (Woosley & Heger 2003) and simulations 
of massive stars exploded by outgoing bi-polar jets (Maeda & Nomoto 2003) show that this 
mechanism may produce substantial amounts of fast moving Ni. It is not clear, though, 
if explosive burning can account for the very fast {v > 0.1c) and massive {M-^i ^ 0.5Mq) 
outflows seen in SN2003dh. Another possibility is that the Ni is synthesized in some sort of 
slow, heavily baryon polluted outflow formed above the black hole (Nagataki et. al 2003). 
The last possibility, investigated here, is that Ni is synthesized in wind-like outflows from 
the accretion disk (MacFadyen & Woosley 1999; MacFadyen 2003). Those authors find, for 
plausible disk viscosities and accretion rates, fast (v ~ 4 x lO^cms"-*^) Ni-rich outflows blown 
off of the inner accretion disk by viscous heating. We provide a simple treatment of these 
winds with an eye toward understanding the wind energetics and the conditions needed for 
efficient Ni synthesis. 

We also discuss the possibility of synthesizing r-process elements in coUapsar events. 
The r-process (Burbidge et al. 1957; Gameron 1957) accounts for roughly half the heavy 
nuclides above the iron group, producing characteristic abundance peaks at A ~ 80, 130, 
and 195. The astrophysical site for the production of these elements remains uncertain. By 
examining the relative isotopic abundances of ^^^Hf/^^'^Hf and ^^^I/^^^I, Wasserburg, Busso, 
& Gallino (1996) argued that at least two distinct r-process sites, operating in the galaxy with 
different rates, must produce the nuclides with A < 130 and A > 130 (see also Wasserburg & 
Qian 2000). Subsequent work by Qian, Vogel, & Wasserburg (1998) and Qian & Wasserburg 
(2000) argues that events which produce the heavy, third-peak nuclei occur ten times more 
frequently than events which produce the lighter (second peak and below) elements. The 
latter must also produce iron-group elements copiously. Gonversely, the high frequency 
(producing A ~ 195) must produce little iron. The association of lower event frequency with 
large Ni production, plus the work of Heger et al. (2003), which argues that the coUapsar 
rate might be ~ 10% of the total core-collapse supernova rate, points potentially to coUapsar 
events as the astrophysical site for ^4 < 130 r-process element synthesis. We argue that the 
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inner coUapsar disk material with low electron fraction may be ejected rapidly, with modest 
entropy, in magnetically dominated filaments or bubbles and that this intermittent outflow 
is a likely site for r-process nucleosynthesis up to >1 ~ 130. 

2. Disk Outflows and Nucleosynthesis 

Nucleosynthesis in outflows from the disk is primarily sensitive to the electron fraction 
(y^), the entropy per baryon (s/kb), and the timescale characterizing the expansion of the 
fluid around the time of efficient a-particlc formation (Hoffman, Woosley, & Qian 1997). To 
estimate these quantities we examine two different, and in some sense limiting, realizations of 
the outflow. The first is a hydrodynamic picture of the outfiow, which assumes the presence 
of stable and ordered pressure profiles: a steady disk wind. The second picture wc examine 
is one where rapid magnetic reconnection or turbulent viscous heating deposits entropy (and 
energy) in localized "bubbles" within the disk, causing rapid ejection of low Y^. material. 
Both types of processes likely occur to some extent. 

3. A Hydrodynamic Wind Picture 

Steady or quasi-steady spherical winds have been extensively used to study nucleosyn- 
thesis in the neutrino-driven wind occurring several seconds after core bounce in core-collapse 
SNe (Duncan, Shapiro, & Wasserman 1986; Woosley et al. 1994; Takahashi, Witti, & Janka 
1994; Qian & Woosley 1996; Cardall & Fuller 1997; Sumiyoshi et al. 2000; Otsuki et al. 2000; 
Wanajo et al. 2001; Thompson, Burrows, & Meyer 2001). Of particular interest to the present 
study is the paper by Qian & Woosley (1996), which provided insight into and analytic ex- 
pressions for describing the connection between the various parameters determining the wind 
itself (neutron star radius and mass, neutrino luminosities, etc.) and the properties of the 
flow that determine the resulting nucleosynthesis: the electron fraction, entropy, and dynam- 
ical timescale. It should be noted that so far no agreed upon neutron star wind solutions give 
a robust r-process (with the possible exception of winds from very relativistic neutron stars 
(Cardall & Puller 1997; Thompson, Burrows, Sz Meyer 2001) or highly magnetic neutron 
stars (Thompson 2003)). 

In this section we describe our simple model of time-independent disk winds, akin to 
the models previously developed in the neutron star context. Roughly, our assumption of a 
steady-state outflow is an attempt to use an estimate of the vertical disk pressure gradient to 
derive the dynamical timescale of the outflow and the increase in entropy due to viscous and 
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neutrino energy deposition. Because we parameterize flow streamlines and currently employ 
only a local a-disk prescription for energy deposition, our analysis is limited and must be 
tested eventually against full MHD simulations of coUapsar disk winds. 



3.1. Equations 

The equations governing the wind structure are determined by mass, momentum, energy, 
and lepton number conservation. For steady flow these are, respectively, 

V-{pu) = (1) 

pu-Vu = f (2) 

u-{V-T) = -u ■ {V ■ T^) + q.pNa (3) 

K-(Vre) = -K-jYe-^-^^)+X..Jl-Y^-^-^^). (4) 



Here p is the mass density, u is the velocity, Qi, is the net neutrino energy deposition rate per 

— * 

baryon, and / is the force (gravitational+viscous+pressure gradient) acting on outflowing 
fluid elements. T — pu<S>u+pg is the stress-energy tensor of the fluid, with p the pressure and 
g the metric tensor. Our treatment is entirely Newtonian, and in spherical polar coordinates 

{1 = r, 2 = 9, 3 = (f)) gu = 1, 5^22 = r^, gas = r"^ sin^ 9, and gij = for i 7^ j. Including 
relativistic effects in aspherical viscous flows is straightforward, but more complicated than 
including relativistic effects for inviscid spherical outflows from neutron stars. 

In eq. (4) we have approximated the flow as consisting of free nucleons and a-particles. 
Ye is the net number of protons per baryon in the flow, Ae-p and \e+n arc the rates for 
electron and positron capture on free nucleons (e.g. Qian & Woosley 1996), and Xfree is the 
mass fraction of free nucleons in the flow. Lepton capture on heavy nuclei is typically very 
slow compared to lepton capture on free nucleons and is neglected here. When Xij-^e = (all 
a-particles) , then, eq. (4) gives zero rate of change for Yg. We approximate Xfj-^e as 

>Tl!lr f -7.074: 



X,„ = 8.2.10»-piexp^^j (5) 

or unity, whichever is smaller (Woosley & Baron 1992). Here Tmcv is the temperature 
in MeV. Eq. (4) does not account for the influence of neutrino capture on Yg. Including 
neutrinos in a reliable way is beyond the scope of this paper. However, in section 3.2.2 we 
outline in general terms how neutrinos are expected to influence Yg. 

The viscous stress tensor is (e.g. Mihalas & Mihalas (1984)) 

Ta = 2/iL>, (6) 
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with 

Here a semi-colon represents the covariant derivative. Following Stone, Pringle & Begelman 
(1999), we neglect all components of D except D'"'^ = D'f^ = (l/2)tij and D^'t' = D'f'^ = 
(l/2r^)w^. The neglect of the other components of the stress tensor is justified if the magnetic 
instabilities providing the shear stresses produce small poloidal stresses. For the coefficient 
of shear viscosity we adopt the parametrization often used in studies of viscous disks: = 
ap/Qk{i^o)- Here a is the standard disk alpha parameter (Shakura & Sunyaev 1973), Tq is 
the radius from which the wind leaves the disk, and the Keplerian frequency is Qk{fo) = 
\/ GM/r^, with M the mass of the central black hole. 

To proceed, we parametrize the trajectories of the outflow by r = (r, 6, 0) = {rQg{6), 6, (p) 
and the wind velocity by -u = (r, 6, (p) = {rog'f, f,^)flk{i^o)- Throughout the paper a prime 
denotes differentiation with respect to 6. We also define the velocity in the r — 9 plane, 
V = f^k{^o)\/ g"^ + g''^. The basic idea of parametrizing the trajectories in the way we have 
comes from the pioneering work of Blandford & Payne (1982). Unhke that work, however, 
we are not solving for the outflow and magnetic field configuration in a global and self 
consistent way. Rather, we will choose trajectories g(d) in order to get estimates of the 
wind parameters important in determining nucleosynthesis. Solving for the trajectories in a 
consistent way likely requires a hydromagnetic simulation which can capture the interplay 
between the outflow and the magnetic fleld conflguration. 

With the above definition for the streamlines, eq. (1) becomes 

pva = Constant. (8) 

Here a = g^sm9{g^ + g'^)~^^^ and is proportional to the area defined by fluid streamlines. 
The equation governing v is found by taking the projection of eq. (2) along the streamlines: 

W 1 - ^ = A - ^ + Pp-, 9 
\ J p a 

where Pp = dP/dp\s — cf, = dP/ds\p, and 

^ = rl^ie sin^ 0^ + rlnl (10) 

is the sum of the gravitational and centrifugal forces. Eq. (3) can be recast as the equation 
governing the entropy of the outflow 
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Here s is the entropy per baryon scaled by kf,. We approximate s as the sum of contributions 
from relativistic hght particles (photons and pairs), free nucleons, and alpha particles: 



s = 5.21^ + Xfree f 12 + In f ^1 I + L^^^ I 15.4 + in f ^ 1 1 , (12) 



where ps — p/lO^gcm ^. The influence of nuclear recombination on the wind dynamics is 
discussed in section 3.2.4. The equation describing the evolution of Ye is 



1 — ^free \ - / - - ^ 1 — Xi 



mr: = -Ae-p ( n - ^^^j + v„ ( 1 - n - ) . (13) 



We adopt a crude parametrization for q^: 

(iv = qh- Qc- (14) 

Neutrino energy loss from the wind occurs principally via capture on free nucleons 
(e~p — s> nz/e and e+n — > pi^e)- The energy loss rate associated with these processes is 
Qc ~ 2. 3rj^gy(MeV/sec baryon). As in the neutrino-driven winds that occur in the late-time 
core-collapse SN cooling epoch, neutrino heating occurs principally via charged- current neu- 
trino capture on free nucleons (z/gn pe^ and Uf.p ~^ ne^). The heating rate for these pro- 
cesses is Qh 5L^,5i((E^)/10MeV)2(l/r^)(MeV/sec baryon) (Qian & Woosley 1996). Here 
we have approximated the inner, neutrino luminous portions of the disk as being spherical, 
Lu,5i is the sum of the i/e and i7e neutrino luminosities in units of lO^^erg/sec, ej, is the average 
i/g or z7g energy, and ry = r/lO^cm. 

To get an idea of the relative importance of viscous and neutrino heating, note that 
for a 3Mq black hole with Kerr parameter a — accreting at O.IMqSCC^^, Popham, 
Woosley, & Fryer (1999) estimate L^^^^i 3 — 7, for a in the range 0.03 — 0.1. Typi- 
cal average neutrino energies for these disks are (E^,) 15 — 20MeV. At r — lO^cm, 
the viscous heating rate is 140T]v[ev(Q;/0.1)(MeV/sec nucleus) for ^ 1, while Qh Ri 
100(L,^,5i/5)(e,y/20MeV)^(MeV/secnucleon). For these disks, then, neutrino and viscous 
heating can be comparable. For a disk surrounding a black hole with high Kerr parameter, 
neutrino heating can dominate over viscous heating. Popham, Woosley, & Fryer (1999) esti- 
mate that a disk with a = 0.95 has a neutrino luminosity about eight times larger than the 
same disk with a = 0. As we discuss though, the parameters important for nucleosynthesis 
are not very dependent on the neutrino luminosity (though the influence of neutrino losses 
is important in determining the disk structure). 
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3.2. Wind Profiles: Yg? s, and Tdyn 

With the above formahsm we can discuss conditions in outflows from the disk and 
imphcations for nucleosynthesis. For simphcity, fluid streamhnes are taken to be straight 
hues making an angle 9o with the plane of the disk (or n — Oq with the z-axis). We assume 
that ^ decreases with distance z above the disk as 

e = exp(-^/e.) (15) 

with the scale height for the decrease in ^. Thus, as the wind moves out of the plane of the 
disk, its velocity in the direction evolves as ^{z)flk- Our discussion and parametrization 
of these steady state winds is similar in some ways to the work of Daigne & Mochkovitch 
(2002), who examined the conditions needed for ultra- relativistic (Lorentz factor much larger 
than unity) outflows from accretion disks. 

We do not present an exhaustive survey of wind models - the simplicity of our model 
probably does not warrant it. Instead we outline how the parameters influencing nucleosyn- 
thesis in the wind depend on the character of the accretion disk and on the starting radius 
of the outflow. In Table 1 we show results from wind solutions for outflows from a moderate 
viscosity {a = 0.1) disk and from a low viscosity (a = 0.03) disk. In both cases the accre- 
tion rate of the disk is M = O.lMosec^^. For each type of disk, outflows from a moderate 
radius {tq = 10^) and a small radius (ro = 10^'^cm) arc considered. Initial (in-disk) values 
of the temperature and density for the calculations were taken from the results presented 
in Popham, Woosley, & Fryer (1999). Typical disk temperatures are a few MeV, typical 
densities are p > lO^gcm"^, and typical entropies are of order 5-10. The starting electron 
fraction was taken from Pruet, Woosley, & Hoffman (2003). The electron fraction in the 
disk depends sensitively on the mass accretion rate and viscosity, and can be anywhere in 
the range 0.1 < l^e ^ 0.53. Results shown in Table 1 were calculated with a neutrino heating 
rate qh = (100/ry)(MeV/sec nucleus). We flnd that changing by a factor of two in either 
direction does not have a big influence on the asymptotic wind parameters. Increasing qh by 
a factor of eight results in an increase in the mass outflow rate of ~ 25% and in increase in 
the asymptotic entropy of 3 units for models C and D. All calculations in Table 1 are for 
9o — 80°. Effects of changing are discussed in §3.3. 

To give a point of reference for the following discussion we show typical wind solutions 
in Figures 1 and 2. Figure 1 corresponds to a wind beginning at ro = lO^cm in a disk 
with a — 0.1 accreting at a rate of O.IMqSBC'^ onto a central black hole of mass 3Mq. 
These parameters are close to those thought to describe conditions in coUapsars. For the 
calculation in Fig. 1, Oq — 80° and = 2ro — 2 x lO^cm. Figure 2 shows a wind for the 
same parameters as in Fig. 1 except with ro = ^2/2 = 10^-^cm. These winds bear qualitative 
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similarities to i/— driven winds from neutron stars. The temperature at the base of the wind 
is approximately that for which the heating and coohng rates balance each other. Also, 
most of the heating occurs at the base of the flow, with the evolution at larger radii being 
isentropic. 



3.2.1. The Asymptotic Entropy 

The trend of greater increase in final entropy with increasing initial gravitational poten- 
tial and the weak dependence of the final entropy on the heating rate (a) is evident. This 
is similar to the case for winds from NS's. Qian & Woosley (1996) argued that the final 

2/3 

entropy should scale as ~ Tq , with only a weak dependence on the heating rate, which is 
a fair approximation to the results shown in Table 1. Overall, the final entropies expected 
for winds from the disk are rather modest (~ 30 — 50) and more typical of the a-process 
than the r-process. 

So far it is not clear to what extent the wind properties relevant for nucleosynthesis are 
determined by our simplified model and to what extent the wind properties are determined 
by more basic considerations. To get insight into the distinction note that the hydrodynamic 
equations can be recast in the form 

Q' = Ts' = h' (16) 

with Q the total energy per baryon added to the flow by viscous and neutrino heating and 
h the BernouUi integral for this system: 

b = -^+TSrad+ ^ + enuc ^(1 + C/c) (17) 

« 500 Q' + TMev.,ad + ^ + enuc -45(1 + f/c)(^)^(MeV). (18) 

Here Srad is the entropy per baryon in radiation, A is the mean atomic mass of bound nuclei 
and enuc is the mean nucleon binding energy. If everything burns to ct-particles the difference 
in enuc before and after nuclear recombination is about 7.074MeV, while if everything burns 
to Fe-group nuclei the difference is about 1.4 MeV larger. The quantity Uc accounts for the 
influence of the (j) component of the velocity in overcoming the gravitational potential. For 
example, within the disk rQ(p ^ rQQk{i^o), so the material already has half the kinetic energy 
needed to escape the gravitational pull of the black hole. In our treatment 

Uc={9/2) f eig'sm'eydO. (19) 

Jn/2 



-9- 



The change in entropy of the outflowing fluid is 

As = Sf-Si^ -, (20) 

where the subscripts i and / denote values at the base and end (i.e. asymptote) of the flow 
respectively If the temperature decreases along outflowing streamlines 

As > Ab/T,. (21) 

Approximate equality in eq. (21) holds if most of the heating occurs very near the base of 
the flow. 

For a given starting radius and disk composition the above considerations give a mini- 
mum value for the final entropy as a function of the asymptotic velocity of the outfiow and 
Uc- Values for the different cases considered in the wind calculations are shown in Table 3. 
The calculations there assume Uc — (or equivalently that is small). In this case the 
centrifugal potential plays no role. There is fairly close agreement between the asymptotic 
entropies found from consideration of the Bernoulli integral and the asymptotic entropies 
found in our wind calculations (Table 1). Entropies found in our calculations are typically 
about 5 units lower than the values in Table 3. This arises mostly because our wind calcu- 
lations have modest and the rotational velocity plays some role in decreasing the effective 
gravitational potential. 



3.2.2. The Asymptotic Electron Fraction 

In contrast to the way in which the final entropy is set, the final electron fraction in 
these disk winds is set by quite different factors than in NS winds. In winds from NS's, 
neutrinos dominate both the energy deposition rates and the lepton capture rates. The neu- 
tron to proton ratio comes into approximate equilibrium with the neutrino spectra. Because 
the neutrinos originate from the neutron rich crust of the neutron star, which has a high 
opacity to electron neutrinos, the Pg spectrum is hotter than the spectrum. Consequently, 
neutrino capture above the nascent neutron star leads to a neutron-rich wind favorable for 
the r— process. 

In winds from accretion disks that are optically thin to neutrinos all factors conspire 
to make > 0.5. In the first place, e"^ capture, rather than neutrino capture, generally 
sets Ye in the disk and in the wind. As viscous heating adds entropy to the outgoing fluid 
the electron degeneracy is removed. Weak equihbrium then favors Ye > 0.5 because of the 
neutron-proton mass difference. 
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Secondly, when neutrino captures are important they tend to increase Ye (Surman & 
McLaughhn 2003; Beloborodov 2003). Very roughly, this can be thought of as a consequence 
of the neutrinos carrying net lepton number away from the deleptonizing disk. To make these 
arguments more quantitative, note that the ratio of the rates for Ug and i?e capture is 



(Qian & Woosley 1996). In eq. (22) A = 1.293MeV is the neutron-proton mass difference, 
{Ely) is the average neutrino energy, and rii, is the neutrino number density. The factor 1.2 
weighting (Ei,)^ is approximate and depends on details of the neutrino spectrum. Taking 
typical neutrino emission as originating from r ~ 10^'^cm in the disk, the disk parameters 
from (Popham, Woosley, & Fryer 1999) and the electron fraction from (Pruet, Woosley, & 
Hoffman 2003), the ratio in eq. (22) can be calculated (Fuller, Fowler, & Newman 1982). 
We find R = 0.94 for the M = 0.1MQSCc-^,a = 0.1 disk and R = 0.66 for the M = 
O.IM0 sec~^, a = 0.03 disk. Because R < 1, neutrino capture above the disk tends to drive 
the outflow proton rich. 

It is worth noting that in the very inner regions of the disk the composition can fall out 
of weak equilibrium (Pruet, Woosley, & Hoffman 2003), with smaller than the equilib- 
rium electron fraction Fe,eg- This results in a low electron Fermi energy, and consequently, 
relatively underluminous Vg emission. For example: for the a = 0.03 disk, material within 
the disk at r = 10^-^ cm has Yg/Ye^eg 0.7 and neutrinos originating from this region are 
characterized by i? = 1.2. By themselves, neutrinos from a region like this would tend to 
drive the outflow neutron rich. However, relativistic effects quash the inffuence of neutrinos 
from the innermost regions of the disk. The neutrino capture rate far above the disk goes 
approximately as the flfth power of the redshift factor h = ^^1 — 2M/rQ (e.g. Pruet, Fuller, 
& Cardall (2001)). This reduces the influence of neutrinos from tq = 10^'^cm by about 
an order of magnitude relative to the influence of neutrinos from r = 10^'^cm. For disks 
accreting more rapidly than pa IMq sec^^, the inner regions become opaque to neutrinos and 
a more careful treatment of the neutrino spectrum is needed (Surman & McLaughhn 2003). 

The above considerations about Ye are exemplifled in Table 1. In all cases the asymptotic 
electron fraction is larger than the in-disk electron fraction. The scaling of flnal Yg with tq, 
^, and a has clear origins. All else being equal, a smaller a implies a denser disk, with faster 
weak interaction rates. In addition, a higher density implies - for a given mass loss rate - a 
lower outflow velocity and more time for weak processes to operate. These are the reasons 
why Yg is so large in the a — 0.03 disk. Similar reasons are behind the scaling of Yg with tq. 
A smaller tq implies a larger density. As well, the entropy of the outflow (and the relative 
importance of positron capture) increases with decreasing tq. Lastly, as increases, the 



R = 



XivgU e-p) 



1.2(£;pJ^ + A^-2A(£;pJ 

i.2(£;,j2 + A2 + 2A(£;,j 



(22) 
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wind material is flung out centrifugally, attaining larger velocities at smaller distances from 
the disk, and there is less time for positron capture. This is why the asymptotic electron 
fraction in the = 4ro case is relatively low. 

Though disks with low a. (or high M) can be very neutron rich (Pruet, Woosley, & 
Hoffman (2003); Beloborodov (2003)), wind-hke outflows will not preserve the neutron 
excess. In particular, wind-like outflows as discussed here cannot result in an asymptotic 
Ye ^ 0.4. An exception to this is for the relativistic jet originating very near the black hole. 
Neutrino-antineutrino annihilation and relativistic effects dominate such ultra-relativistic 
outflows and they can remain neutron rich (Pruet, Fuller, & Cardall 2001). 

It should be noted that we may overestimate the electron-fraction in the a = 0.1 disk. 
This is because our calculation shows about half of the change in Y^. coming within one- 
pressure scale height of the disk mid-plane. Though our calculations show the outflow to 
be in near hydrostatic equilibrium near the disk, a two-dimensional calculation of the disk 
structure and composition would give a clearer picture of how evolves in outflows from the 
disk. The uncertainty in Y^ is unfortunate because efficient ^^Ni synthesis hinges sensitively 
on Yf. being larger than 0.5. For the a — 0.03 disk, the asymptotic Y^ is less sensitive to the 
vertical disk structure (at least in our simple calculations) because most of the change in Yf. 
occurs a few pressure scale-heights above the disk mid-plane. 

3.2.3. The Dynamical Timescale 

Both the final entropy and electron fraction are set by processes near the disk. By 
contrast, the dynamic timescale at the epoch of nucleosynthesis is determined by the wind 
structure at r ~ 2 — 5 x lO^cm. One way to estimate the timescale characterizing the 
expansion of the fluid at T < l/2MeV is simply to use the calculated wind proflles. This 
is likely not correct. It seems implausible to expect that the disk outflow will remain well 
collimatcd in quasi-cylindrical geometry for z > lOrg. More likely is that the magnetic or 
pressure conflnement breaks down at large radii and the wind assumes a quasi-spherical 
expansion and begins to coast. A further dynamical affect may also influence rdyn at large 
radius: when outflow encounters the overlaying stellar mantle, which also happens in NS 
winds, the wind will be slowed. 

For lack of a calculation of the interaction between the disk wind and the exploding 
star, we make the assumption that after the sonic point the wind begins expansion with 
a oc (here a is the area deflned in eq. 8). While the wind has appreciable enthalpy the 
expansion is roughly homologous with v ca r and a dynamic timescale Tdyn — r/v. Entropy 
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and mass conservation imply the scaling T oc p^/^ oc 1/r. Once the velocity asymptotes to 
V{ a coasting phase described by pvfr^ = const, follows. In the coasting phase the dynamic 
timescale is again approximately r/wf, though it is now an increasing function of radius. In 
Table 2 we list the dynamic timcscales calculated as described here. The first Tdynamic listed 
for each calculation is that appropriate for homologous expansion (rsomc/'^somc)- The second 
dynamic timescale listed is that for the flow when it is coasting. For consistency with the first 
definition of Tdyn this dynamic timescale is defined as the time needed for the temperature 
to decrease by a factor of e from Tg = T/IO^K — 5. The two different dynamic timescales 
should approximately bracket the plausible range of dynamic timescales. 

Determining the dynamic timescale is equivalent to determining the mass outflow rate 

M(sph) ^ 4T^r^pv^ (23) 

for a given asymptotic entropy. Here M^^p*^) is the mass outflow rate that would obtain if 
the outflow were spherical. For winds from the inner regions of accretion disks the true mass 
outflow rate is typically much smaller that M^'^p'^) because of the coUimation of the wind. 
Observations of ^^Ni from GRBs may help constrain M'^'^p'^^ (see below) . Taking r ~ Vft gives 
the time at which a given temperature is reached in the outflow 



t = 0.022 sec \l ( -^f^] ■ (24) 



'^O.l 



Here M_i = M^'p^VO-I^© sec'S S30 = s^ad/30, andvo.i = v/O.lc. Equation (24) determines 
the coasting dynamic timescale (i.e. the temperature e-folding time when Tg = 5) as Tdyn = 



O.lsec yM_iS3o/wg.i. 

At some point the wind from the disk will encounter the overlaying stellar mantle. This 



will influence nucleosynthesis in the wind if the wind is slowed beforer ^ 2xl0®cmy M-iSao/vo.i, 
where Tg ^ 2.5. At times less than a few seconds after disk/black hole formation, then, a 
simple wind picture is inadequate. For winds leaving the disk at times greater than a few 
seconds after disk formation the wind/envelope interaction is not important if the wind is 
energetic. Calculations by MacFadyen & Woosley (1999) show that ait — 9.48 sec the wind 
has cleared a region out to r > 3 • lO^cm in the star. Calculations by Maeda & Nomoto 
(2003) of outgoing bi-polar jets - whose influence on the star may roughly approximate the 
influence of the disk wind - show that by t — 5 sec outflow along the rotation axis continues 
uninhibited out to r ~ 2.5 • lO^^cm. This is somewhat faster than the shock velocities typical 
of ordinary core collapse SNe (vghock ~ lO^cms"^). 

Though the wind/envelope interaction is not expected to influence nucleosynthesis in 
the wind, interaction of the wind with the stellar mantle will ultimately slow the outflow 
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and influence the observed Ni velocity. As a very rough estimate, if 2Mq of Ni-rich ejecta 
mixes with IOMq of stellar mantle, the observed Ni velocity will be a factor of •\/2/10 ~ 1/2 
smaller than the estimates in Table 1. 



3.2.4- The influence of nuclear recombination on the wind dynamics 

During a recombination the total entropy is constant (apart from the influence of ex- 
ternal heating sources), but entropy is transferred from the nucleons to the e^/7 plasma. 
The amount of entropy transferred can be seen from cq. (12). For the modest entropy 
outflows discussed here. ln(T]^j(,y/p8) ^ 3 during recombination. So, the pair plasma gains 
about (12 + 3) — (1/4) (15.4 + 3) ~ 10 units of entropy. An equivalent way to estimate 
the entropy transfer is to note that a— particles arc bound by ~ 7.074MeV/nucleon, so 
As 7.07MeV/Trec ~ 10, where Tree ~ 0.7MeV is the recombination temperature. ^^Ni 
is bound by 8.6 MeV/nucleon. This means that synthesis beyond He will add more energy 
to the pair plasma. Our calculations neglect the influence of this extra energy input on 
the wind dynamics. In part, this is because of the difficulty of coupling a nuclear network 
to our calculations, and in part because the lion's share of the energy release is from a- 
recombination (7 out of 8.6 MeV). 

Roughly speaking, nuclear recombination can influence the wind dynamics in two dif- 
ferent ways. If recombination occurs below the sonic point there is the potential for the shift 
in entropy (as well as pressure and energy) from nucleons to the pair plasma to change the 
amount of viscous and neutrino heating suffered by the outgoing wind. This would change 
the asymptotic entropy, the mass outflow rate, and the electron fraction in the wind. Of 
the models we discuss in this paper, only the a = 0.03 disks have outflows with sonic point 
temperatures below the recombination temperature. For these we have calculated wind so- 
lutions with and without alpha recombination. The winds without alpha recombination are 
equivalent to assuming that alpha particles are unbound. As can be seen from Tables 1 and 
4, alpha recombination has a modest influence on our estimates of the asymptotic entropy 
and mass outflow rate. This is not so surprising, since recombination occurs at relatively 
low temperatures and after most of the heating has occured. 

The second, and most important effect of recombination is simply to make the nuclear 
binding energy available for kinetic energy (eq. 17). Indeed, as Woosley & Heger (2003) 

have pointed out, the observed kinetic energy for SN2003dh (~ 2.6 — 4 • lO^^erg) is not so 
different from the energy liberated by recombination of a solar mass of free nucleons to Ni 
(~ 1.65- lO^^erg). Observationally, though, it still cannot be determined if Ni recombination 
accounts for essentially all of the observed kinetic energy, or only ~ 1/3 of the observed 
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kinetic energy. Our simple calculations indicate that Ni recombination accounts for about 
70% of the kinetic energy of SN2003dh if the accretion disk has a. — 0.03, and about 25% of 
the kinetic energy if the accretion disk has a = 0.1. 

3.3. Outflows Flung Magnetically from the Disk 

In the above discussion we assumed that the bulk of the work in driving the outflow 
is done by pressure/entropy gradients established by viscous and neutrino heating. It is 
also possible that material can be flung outward along a magnetic field line with little or 
no help from pressure gradients. This case has been thoroughly discussed by Blandford & 
Payne (1982), who showed that such outflows might mediate angular momentum transfer 
in tenuous accretion disks. In collapsar environments it is unclear to what extent such a 
mechanism can operate because the outflows arc so dense and their inertia is important. 

It is worth noting, though, that outflows centrifugally pushed along a magnetic field line 
can be qualitatively different from the wind-like outflows discussed above. In general they 
are neutron-rich and do not synthesize ^^Ni. The asymptotic entropy and electron fraction 
can both be very low - even lower than the in-disk values if neutrino cooling is important. 
For example, if for model D is changed to 70°, then the asymptotic entropy and Y^. become 
23 and 0.46 respectively. This can be understood from the discussion of the Bernoulli 
integral in section 3.2.1. For = 70° and — 2ro, the rotational velocity at 2; = 2ro is 
about 50% larger than the kcplerian velocity there, so that thermal heating does not have 
to do much work (i.e. 1 -|- t/g is small). If = 70° and = l^'O) so that heating must do the 
work against gravity, the asymptotic entropy and electron fraction are essentially the same 
as for = 80° (s = 28, Yg = 0.50). To summarize, the larger that magnetic or other effects 
keep r(j), the smaller will be the final entropy and electron fraction. 

3.4. Nucleosynthesis in Winds 

Here we concentrate on the synthesis of radioactive ^^Ni. There is no significant pro- 
duction of r— process elements for the high electron fractions and modest entropies found in 
our calculations. Although material leaving the disk from very near the hole will have higher 
entropy, Yg will be too large for the r— process except perhaps in outfiows with very rapid 
expansions (Meyer 2002). 

Winds that are accelerating at the epoch of nucleosynthesis (T ~ 0.5 MeV) generally 
expand too quickly for efficient Ni production. There is no time for the 3-body reactions 
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that lead to ^^C synthesis and efficient a— captures. In contrast, the evolution and nucle- 
osynthetic yields of winds that are not accelerating, but coasting, are largely described by 
two parameters (see eq. 24). For convenience we take these to be the asymptotic entropy s 
and 

P^M_1^ rl^, (25) 

■^0.1 

Ni synthesis also depends on Y^, though the final Ni mass fractions do not vary greatly for 
0.5 < < 0.53. 

In Figure 3 we show final Ni mass fractions as a function of the entropy and the param- 
eter (3. Large Ni mass fractions are favored by lower entropies and larger /3's. At an entropy 
of 50, which is obtained in material leaving the disk from Tq < 3 x lO^cm, Ni synthesis is 
inefficient unless M_i/vq i > 4. For material with an entropy of 30, which is characteristic 
of material leaving the disk at r ~ lO^cm, Ni synthesis is efficient {X^ > 0.25) as long as 
M-i/vl^ > 1/4. Wc note again that the M here is the isotropic equivclant mass outflow 
rate. The true outflow rate is smaller because the Ni-wind is confined to a fraction of the 
solid angle above the disk. 

As far as Ni synthesis in coUapsars is concerned, our results are promising. Table 4 lists 
the mass outfiow rates and final Ni mass fractions for the different wind calculations. Note 
that the mass outflow rates in Table 4 are representative of the "true" mass outflow rates, and 
that M'^^P'^) (deflned through eq. 23) is much larger. The outflows we calculate for material 
leaving the disk at lO^cm all efficiently synthesize Ni and have rather fast expansion velocities 
V > 0.1c. However, the simple considerations presented here suggest that a low-alpha disk 
may have trouble synthesizing ~ O.5M0 of Ni. The reason is that such a disk is efficiently 
cooled by neutrino losses, so that the disk material is tightly gravitationally bound (see the 
BernouUi parameters in Table 3), and heating cannot drive large mass outfiow rates. For the 
a — 0.03 disk, the mass outfiow rates are about an order of magnitude smaller than for the 
a — 0.1 disk. For the a = 0.1 disk our models predict outflow rates inr^poVo ~ O.O7M0 sec~^ 
for tq = 10^ and a factor of about 3 smaller for material leaving the disk from tq = 10^'^cm. 
These numbers are in the right range for explaining a total disk cjccta mass of ~ 2Mq from 
an event with a duration typical of long duration GRBs (~ 10 — 100 sec). 



4. Impulsive Rapid Mass Ejection: Magnetic Bubbles 

We imagine a background disk wind, similar to the solutions obtained in the preceding 
section, but modulated by rapid impulsive events due to magnetic reconnection. Highly 
magnetic (low rj = P/(5^/87r)) filaments may be formed in the disk mid-plane and emerge 
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into the steady wind background rapidly as they expand and accelerate in the approximately 
exponential atmosphere of the disk. We propose here that this environment will have modest 
entropy, but low electron fraction (similar to the that obtains at the disk mid-plane) on 
account of the rapid expansion of the bubbles. It is in these events that A ~ 130 r-process 
elements may be synthesized. 

There is an important difference between heating from magnetic instabilities and heating 
in an alpha-disk model. The viscous heating rate in an alpha-disk tracks the density. In 
eq. (11), for example, it is seen that s ~ P/{pT) ~ Constant, so that very rapid increases in 
entropy can not be obtained in such a model. By contrast, heating from magnetic instabilities 
depends on the geometry of the magnetic fields. There is nothing to prevent heating in a 
relatively baryon dilute region (say two pressure scale heights above the disk mid-plane). 
Such heating can result in rapid, large increases in the entropy of the outflowing material. 
This basic idea is behind a number of suggestions for the origin of relativistic outflow in 
GRBs (Narayan, Paczynski, & Piran 1992; Kluzniak & Ruderman 1998). 

To see how bubbles might form in the coUapsar accretion disk, consider the following 
argument adapted from Kluzniak & Ruderman (1998), who studied differentially rotating 
neutron stars. For simplicity, consider an isothermal disk with a structure and equation of 
state given by 

p{z) = poexp(-2;) (26) 
s{z) = T=^/Vp = soexp(z) (27) 
P{z) = pT = Poexp{-z). (28) 

Here we have assumed the gas to be dominated by free nucleons in the disk and have scaled 
the height above the disk in units of the disk scale height H. A fluid clement rising adia- 
batically from the mid-plane of the disk will be in pressure equilibrium with the background 
fluid. If there is no magnetic pressure, then this element (denoted with the subscript b) will 
have a density 

Pb = poexp(-3z/5) > p{z), (29) 

which implies that the fluid element will fall back to the mid-plane. If the fluid element 
carries a magnetic pressure Pm, pressure equilibrium with the background disk implies 

p \ 3/5 



Pb = Po exp (-3^/5) - J . (30) 

If the magnetic pressure is large enough then pf, < p{z) and the fluid element will be buoyant. 
The force per unit volume on the fluid element is ^\z, which gives the equation of motion 



(31) 

Pb 
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Here we have scaled time in units of Jl^.^. If we make the simple assumption that the 
magnetic pressure increases linearly with time and Pm — vPoi't/'^'^), then eq. (31) becomes 



z — z 



exp (-3^/5) (^1 - ^ exp(^)^ 



-3/5 



(32) 



This equation can be solved for z{t) to provide an estimate of Pm{z) and the entropy increase 
of the bubble as a function of the reconnection height. As a simple approach, suppose that 
an element rises one pressure scale height per radian that the disk rotates, or z — t va. our 
notation. This is a reasonable assumption since the magnetic pressure quickly becomes large 
compared to the disk pressure as the filament rises, li z ^ the ratio of magnetic energy 
density to thermal energy density evolves approximately as 

^ = ■ (33) 

PfeTb 27rexp(— — 772; 

If the energy in magnetic fields is transferred to thermal energy of the buoyant bubble the 
fractional increase in entropy of the bubble is approximately the ratio given in eq. (33). As 
an example, suppose that the initial magnetic field is 10% of the equipartion field (77 = 0.1). 
If the energy in magnetic fields is transferred to the bubble at 1 disk-scale height above the 
mid-plane the fractional increase in entropy will be ^ 6%. If the energy is transferred at 
three scale heights above the mid-plane, entropy will increase by a factor of ~ 23. 

To assess if the above argument seems plausible, let us assume that the magneto- 
rotational instability (the MRI; Balbus & Hawley 1994; Balbus & Hawley 1998) operates in 
collapsar disks and that this produces the turbulent viscosity necessary for accretion. Then, 
the local saturation magnetic field at a radius Tq will be in rough equipartition with the 
azimuthal kinetic energy density; 

^-i-J (34) 

where = roO^ = ro(GM/rQ)^/^. Empirically, for simulations of the MRI in accretion disks 
(Hawley, Gammie, & Balbus 1996) the magnetic field saturates at a sub-equipartition value, 
typically l/27r times the field estimated from eq. (34). Taking this into account, 

S,at(ro) ~ 1 X 10^4 G pI^' mI'^ r-^'\ (35) 

where pg = p/10^ g cm~^, r^^ = tq/IO'^ cm and M3 = M/3Mq. This magnetic field strength 
is in the right range for explaining modest entropy increases in the disk material. If the 
MRI is operating, the time to build a magnetic field is set by f^^^. The timescale for the 
maximum growing mode is 

dQk 



TMax = 47r 



dlnr 



(36) 
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Roughly, we may assume that the local magnetic field may be entirely built and dissipated in 
TMax- Using eq. (34), and estimating the reconnection time as trcc ^ ^Max — L^iec/vA, where 
v\ is the Alfven speed, we find that the characteristic length scale for reconnections Ln^c ^ 
115 km. This may seem uncomfortably large, but we note that the trcc above is the minimum 
required - we should perhaps require 3-4 doublings of the field at ro before dissipating. In this 
case, the reconnection time would be longer and L Rec would be commensurately smaller. We 
expect these reconnection events to happen impulsively, and be separated in time by ~ TRec- 

A bubble will begin to rise after entropy is deposited in a localized region. If the bubble 
is formed in a region of the disk where approximate hydrostatic equilibrium obtains, the 
bubble height z will evolve as 

Z Ri ZoQkt- (37) 

Here wc have made use of the fact that the gravitational force in the vertical direction is 
approximately fllz. The characteristic length scale zq is taken here to be the pressure scale 
height H in the disk. The rising bubble will be in pressure equilibrium with the background 
fluid. This implies that initially the temperature evolves according to 

T oc exp(-^/4ii') (38) 

provided that the bubble is radiation dominated. These considerations suggest that the 
timescale characterizing the initial expansion of the bubble should be of order 

'^initial 

lO/Qk ~ 0.01 sec. Noting that the positron capture rate is Ae+n ~ llsec~^(T/3MeV)^ gives 
an estimate of the change in of the outgoing bubble 

AYe ^ 0.2 ( ' ( ""^"'^"^ ) . (39) 
\3MeY J VO.Olsec/ ^ ^ 

Here Tq is the temperature of the bubble formed after magnetic reconnection. We note 
that for the a = 0.03, M = O.lAfosec^^ disk, the height averaged disk temperature is 
~ 2.2 — 2.8MeV for ro < lO^cm. Bubbles can retain low values of Ye. 

Equations (37) and (38) are appropriate for characterizing the bubble passage through 
the region of the disk that is in approximate hydrostatic equilibrium. Those equations are not 
appropriate for estimating the expansion timescale at the late times and low temperatures 
important for nucleosynthesis. The expansion timescale at T ~ O.SMeV depends on the 
proflle of the background wind. If the background wind is radiation dominated well above the 
disk, then pressure equilibrium between the bubble and wind implies that the temperatures 
in the wind and bubble are equal. In this case the expansion timescale of the wind (r in 
Table 2) may be used as a rough indication of the expansion timescale of the bubble. 



-19- 



4.1. Nucleosynthesis in Bubbles 

As discussed in the introduction, there are indications that the site of the A < 130 
r-process nuclei is different from the site of the r-process nuclei with A > 130. In particular, 
observations of ultra-metal poor stars (e.g. Sneden et al. 1996; Burris et al. 2000) and 
inferences from elemental abundances in pre-solar meteorites (Wasserburg, Busso, & Gallino 
1996) suggest that the ^^^'^^'''Ba isotopes are not significantly produced in the same events 
that produce i27,i29j g^j^^j lighter r-process elements. An exhaustive survey of how a nuclear 
burning site can produce a near solar abundance pattern of nuclei near one of the r-process 
peaks, while leaving the other peaks unpopulated, is beyond the scope of this paper. Here 
we show in broad stroke that conditions in coUapsars are favorable for the synthesis of the 
A ~ 130 r-process elements. In the spirit of investigation, we also discuss a way in which a 
near-solar abundance pattern of (only) A < 130 r-process elements can be synthesized. 

Neutron-rich outflows characterized by rapid expansion generically synthesize r-process 
elements (e.g. Hoffman, Woosley, & Qian (1997)). As an example, we show in Fig. 4 
results of nucleosynthesis calculations for outflows with Yg = 0.2, a temperature e-folding 
time r = 0.12 sec, and three different entropies. These conditions are close to those we 
estimate may obtain in outflows from the inner regions of accretion disks with a — 0.03 and 
M — O.lM0sec~^. Here the (unnormahzed) overproduction factor for nucleus j is defined 
as 

OU) = (40) 

where Xj is the mass fraction of the nucleus j in the bubble and Xqj is the mass fraction 
of the nucleus in the sun. Though it is evident from Fig. 4 that some r-process elements are 
synthesized, good agreement with the solar abundance pattern is not obtained. 

It is unreasonable to expect that a single type of bubble will prevail in the dynamic 
coUapsar environment. Rather, a broad spectrum of bubbles, with different entropies, dy- 
namic timescales, electron fractions, and so on will be created. To investigate the average 
nucleosynthesis we generated random bubbles with properties defined according to 



s = 50 + 50r (41) 

Ye = 0.15 + 0.25r (42) 

T = 0.03(1 + 4r) sec (43) 

Tg,^i. = l + 2r. (44) 



Here r is a random number between and 1 and is generated separately for each of eqs. (41- 
44). The choices for Ye and r above were adopted because they represent the expected range 
of conditions in high entropy outfiows from the inner regions (r < lO^cm) of a disk with 
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a — 0.03 and M — O.lMesec"^. Choices for s and Tg^mix are not well constrained. The 
parameter Tg mix is the temperature at which the bubble is assumed to mix with the proton- 
rich ambient medium. When the bubble shears or destabilizes, the free neutrons available 
for capture are diluted. To represent this wc made the rough approximation that all neutron 
captures then cease and the abundance pattern of neutron rich elements is frozen in (except 
for /3 decays) for Tg < Tg^mix- Because charged particle reactions are relatively slow for Tg < 3 
we expect this to be a fair approximation. 

Average overproduction factors for one hundred bubbles generated as described above 
are shown in Fig. 5. Overall there is quite good agreement with the solar abundance pattern 
of 90 < A < 130 r-process elements. Overproduction factors for ^"°Mo/°'' Ag/^° Pd/^^ Cd/^^ Sb 
and ^^^Sn - elements principally or entirely synthesized in the r-process - are all within a fac- 
tor of a few of the overproduction factor for ^'^^Pd. ^^^I is underproduced by a factor of about 
5. Most of the A ~ 130 nuclei have progenitors about 8 units from stability near the N=82 
closed neutron shell in our calculations. Uncertainties in the location of the shell closures 
and the (5 decay rates of near-drip line nuclei may account for the modest underproduction 
of ^^''l. ^^^Ba and heavier elements are absent. In Fig. 6 we investigate the sensitivity of nu- 
cleosynthesis to the distribution of s and Tg min in the bubbles. Nucleosynthesis is relatively 
insensitive to modest changes in cqs. (41-44). 

Our calculations indicate that typical overproduction factors for 2nd-pcak nuclei are 
O ~ 5 ■ 10'^. With this, we can estimate how much mass collapsars have to eject as bubbles 
in order to account for the galactic inventory of A < 130 r-process nuclei. For events with 
the frequency of Type II SNe, studies of galactic chemical evolution indicate that a typical 
nuclide must have an overproduction factor 

^ (45) 

in order to account for the observed solar abundance of that nuclide (e.g. Mathews, Bazan, 
& Cowan (1992)). In eq. (45) M^-^' 10 - 2OM0 is the typical mass ejected by a type II 
SN and M^X{j) is the mass of the nuclide j ejected. For collapsars, 0{j) must be ~ l//c 
larger, where fc is the fraction of core collapse SNe that become collapsars. The total mass 
ejected as bubbles must then 

If typical bubbles are formed with initial radius r^, temperature T^, and entropy per 
baryon, eq. (46) implies that the number of bubbles needed per event is 

Sb fW^cmV f2McYy Mb 
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Now, if the disk lasts for a time t, the number of disk revolutions, or magnetic field windings, 
per bubble is 

^^'^ 27rnb 50 sec lO'^scc^^ ' 
which is a reasonable number if magnetic instabilities take a few rotations to develop. 

We have argued that if collapsars are only 1% as frequent as Type 11 SNe then each 
coUapsar needs to eject IO^^Mq of r-process bubbles in order to account for the total mass 
of 130-peak r-process material present in the galaxy today. This required ejecta mass is 
approximately 0.01% of the total mass ejected in a Ni wind and is roughly consistent with 
a description of bubbles as forming on a magnetic instability timescale and observed GRB 
durations. However, there still remains the more delicate issue of whether collapsars as the 
130-peak site is consistent with observations of metal poor stars and what is known about 
galactic chemical evolution. 

Unfortunately, much less is known about the production of the A<130 elements than is 
known about the production of the heavier r-process elements. For example, while CS22892- 
052 and a handful of similar stars are thought to have been enriched by only one A > 130 
production event, these authors are not aware of any stars thought to have been enriched by 
only a single A < 130 enriching event. Within the context of the models developed by Qian 
and Wasserburg this has an explanation in the relative rarity of the events which produce 
the second peak elements. 

One type of argument which has been used to differentiate proposed r-process sites is 
based on an analysis of the refreshment rate of material in the interstellar medium (ISM). 
This argument has previously been used to assesss the viability of neutron star mergers as 
an r-process site and can also be applied to collapsars. Qian (2000) showed that enrichment 
of the ISM with second peak elements every lO^yr is consistent with the observed trend of 
Ag abundances with [Fe/H] in metal poor stars. Here we simply note that a refreshment fre- 
quency of /ref ~ (lO^yr)"^ is consistent with enrichment of the ISM by collapsars. Following 
the notation of Qian (2000), 

^'•'-('°V)-(pg^)(4)(^). (49) 

where f^^ is the rate of type 11 SNe in the galaxy and Afmix is the total mass swept up by 
a collapsar remnant. Note that M^i^ for collapsars is about an order of magnitude larger 
than the mass swept up by type 11 SNe. This is because the kinetic energy in the collapsar 
explosion is about ten times larger and because the swept up mass is proportional to -E^/Jetic 
(Thornton et al. 1998). 
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5. Results and Conclusions 

We have considered the nucleosynthesis which may attend outflows in the collapsar 
Gamma Ray Burst environment. Both wind-like outflows and bubble-like outflows were 
considered. Wind-like outflows may be relevant for recent observations of SN1998bw and 
SN2003dh, which hint at a robust connection between the central engines of GRBs and 
core-collapse SNe - or at least SN-hke light curves. To power such hght curves requires 
~ 0.5Mq of radioactive Ni moving outward very rapidly {v > 0.1c for SN2003dh). These 
are characteristics beyond the reach of canonical SNe with energies ~ lO^^erg. The results 
of our simple models of viscosity and neutrino-driven winds are promising. Under a broad 
range of conditions such winds copiously produce fast moving radioactive Ni. 

In general, winds from collapsar disks cannot preserve a large neutron excess. This 
implies that these winds will not synthesize interesting neutron rich elements. However, 
chaotic heating or buoyant magnetic fllaments in localized regions in the disk result in bubbles 
which rise on a timescale comparable to a Kepler period. This is fast enough to preserve the 
neutron excess found in the mid-plane of the disk. Though we do not have a complete theory 
of bubble production, we have shown that the solar abundance pattern for 90 < A < 130 r- 
process elements can naturally be produced in bubble- like outflows and that the requirements 
for the total mass ejected are plausible. As we have discussed, identification of coUapsars as 
the source of the 2nd peak r-process elements is consistent with a number of observational 
indications. 
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Fig. 1. — Wind for material beginning at Tq = 10 cm from a disk with a = 0.1 and 
M — O.lM0sec~^. The sohd hnes are for calculations including a— recombination, while 
the dashed lines are for calculations which neglect the influence of a recombination on the 
wind. Note that in our simple calculations most of the heating (entropy change) occurs near 
the base of the flow. 
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Fig. 2. — The same as Fig. 1 except for material beginning at Tq = 10^-^cm. Note that this 
figure has logarithmic spacing on the vertical axis. 
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Fig. 3. — Plot of Ni mass fraction vs. the parameter (3 for different entropies. This plot 
assumes Yg = 0.51. 
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Fig. 4. — Overproduction factors for nuclei synthesized in bubbles with Yg = 0.2, r = 0.12 sec, 
and different entropies. Solid lines connect isotopes of a given element. The most abundant 
isotope in the Sun for a given element is plotted as an asterisk. A diamond around a data 
point indicates the production of that isotope as a radioactive progenitor. Though some 
r-process elements are synthesized, agreement with the solar abundance pattern is poor. 
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Fig. 5. — Average overproduction factors for 100 bubbles generated according to eqs. (41-44). 
Agreement with the solar abundance pattern of r- process elements with A < 130 is quite 
good, though ^^^I is underproduced by a factor of about 4. Production of species heavier 
than A — 130 is negligible. 
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Fig. 6. — Influence of changing assumptions about the distribution of bubbles produced. 
Run (a) is the reference run and corresponds to the distribution of bubbles described by 
eq. (41-44). Run (b) is the same as run (a) except the entropy in the bubbles is assumed to 
be uniformly spread from s = 50 to s = 150. Run (c) is the same as run (a), except Tgmin is 
assumed to be uniformly spread from 1.5 to 3.5. In run (a) the average of 100 bubbles was 
taken while in runs (b) and (c) the average of 30 bubbles was taken. 
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Table 1. Wind characteristics for different parameters 



Model 


a 


ro(cm) 




Y 

e,i 








A 


0.03 




2 


0.12 


0.56 


54 


0.21(0.17] 


B 


0.03 


10^ 


2 


0.22 


0.50^ 


24^ 


0.11[0.06] 


C 


0.1 


106.5 


2 


0.44 


0.54 


52 


0.36[0.34] 


D 


0.1 


10^ 


2 


0.43 


0.50 


33 


0.23[0.19] 


E 


0.1 


106.5 


4 


0.44 


0.50 


47 


0.32[0.30] 


F 


0.1 


10^ 


4 


0.43 


0.46 


32 


0.22[0.18] 



^For model B sy^ = 26 and Yg,/ — 0.53 if alpha recombination 
is neglected. For all other models alpha recombination does not 
effect s/ or Y^j in our calculations. 

^Values in square brackets are those calculated if alpha recom- 
bination is neglected. 
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Table 2. Dynamic timescales for the different winds. 



Model 


^homologous (sec) 


''"coast (sec) 


A 


0.03 


(NC)^ 


B 


0.17 


(NC)^ 


C 


0.007 


~ 0.04 


D 


0.03 


~ 0.10 


E 


0.002 


~ 0.03 


F 


0.012 


- 0.08 



*The sonic-point temperature for 
these winds is less than O.SMeV. 
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Table 3. Minimum values of the final entropy for a given final velocity 



a 


ro(cm) 




S/l'yfinal = O.lc) 


S/(^^final = 0.2c) 


S/(^^final = 0.3c) 


0.03 




-132 


53 


58 


66 


0.03 




-37 


26 


33 


44 


0.1 




-83 


45 


49 


56 


0.1 


10^ 


-21 


28 


35 


46 



^In-disk value of the BernouUi parameter (eq. 17). 
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Table 4. Mass ablation rate and estimated final ^^Ni mass fractions for some different 

wind calculations. 



Model 


M{Mq sec-^)^ 




A 


10-=^ 


<o.i 


B 


4 X 10-^ 


~ 0.5 


C 


3 X 10-2 


~ 0.1 


D 


7 X 10-2 


~ 0.4 


E 


6 X 10-2 


<o.i 


F 


0.14 


_c 



^Defined as Anr^poVo. 

•^Estimated from Figure 3 under 
the assumption that the wind be- 
gins quasi-spherical expansion after 
the sonic point. Dynamic timescales 
were taken from Table 2. 

^This model has < 0.48 and 
does not produce ^^Ni. 



